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Quantum Monte Carlo method is used to study the coupled spin-pseudospin Hamiltonian 
in one-dimension (ID) that models the charge-ordering instability of the anisotropic Hubbard 
ladder at quarter filling. We calculate the temperature dependence of the uniform spin sus- 
ceptibility and the spin and charge excitation spectra of the system to show that there is a 
parameter and temperature region where the spin degrees of freedom behave like a ID anti- 
ferromagnetic Heisenberg model. Anomalous spin dynamics in the disorder phase of a typical 
charge-ordered material a'-NaV20s is thereby considered. 
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1. Introduction 

Charge-ordering (CO) instability has recently been 
one of the major topics in the field of strongly corre- 
lated electron systems. Here, elucidation of the observed 
anomalous behaviors of electrons associated with the CO 
phase transition has been the central issue. This includes 
questions on the charge dynamics above the transition 
temperature Tco as well as on the CO spatial patterns 
realized below Tco- A well-known example is the vana- 
date bronze a'-NaV^Os where the system may be mod- 
eled as a lattice of coupled ladders (or a trellis lattice) at 
quarter filling. 1-5 ) Strong intersite Coulomb interaction 
between electrons is believed to be the origin of the CO 
instability. 2 ' 3 ^ In this material, the CO with a zigzag or- 
dering pattern is observed below T co = 34 K, 6-10 ) and 
associated with this, a number of anomalous behaviors, 
which can be related to the slow dynamics of charge car- 
riers (or charge fluctuation), have been observed above 
7co- 10 ~ 17) Anomalous response of the spin degrees of 
freedom has also been noticed. 9, 18, 19 - ) It seems therefore 
quite natural to wonder how in such systems the spin 
degrees of freedom behave near the CO phase transition 
when they are on the slowly fluctuating charge carriers. 
In this paper, we thus consider the issue: what are the 
consequences of charge fluctuation at T > Tco to the 
spin degrees of freedom? 

One of the simplest models that allow for such situa- 
tion is the anisotropic Hubbard ladders at quarter filling 
with the strong intersite Coulomb repulsion. We here use 
an effective Hamiltonian written in terms of the spin and 
pseudospin (representing charge degrees of freedom) op- 
erators. 5, 17 > 20 > 21 ) This Hamiltonian is derived from the 
Hubbard ladder model by the perturbation theory 5, 20, 21 ) 
where the hopping parameter between the rungs of the 
ladder is assumed to be small compared with the onsite 
and intersite Coulomb repulsions as well as the hopping 
parameter in the rung (i.e., the anisotropic ladder). 3 ) Al- 
though the CO is not realized in this model (since it is 
the ID quantum-spin model), we can simulate anoma- 



lous behaviors of the spin degrees of freedom under the 
strong charge fluctuation. We will apply the quantum 
Monte Carlo (QMC) method to this model to calculate 
the temperature dependence of the uniform spin sus- 
ceptibility and the spin and charge excitation spectra, 
thereby clarifying consequences of the interplay between 
its spin and charge degrees of freedom. 

We note that, in this coupled spin-pseudospin model, 
the spin exchange interaction is necessarily associated 
with the charge excitation; i.e., the spin excitations can- 
not occur without making the exchange of the pseu- 
dospins. We will then show that nevertheless there 
is a parameter and temperature region where the spin 
degrees of freedom behave like a ID antiferromagnetic 
Heisenberg model; i.e., the spin degrees of freedom are 
'separated' from the charge degrees of freedom in this 
region. We will moreover show that the spin system 
behaves in different manner depending on whether the 
temperature T is below or above a crossover tempera- 
ture T* which is related to the pseudospin excitations; 
at T < T*, it behaves like a ID antiferromagnetic Heisen- 
berg model with a T-independent effective exchange cou- 
pling constant J e ff with large renormalization, whereas at 
T > T*, J c ff decreases rapidly with increasing T, where 
the effective Heisenberg-model description ceases to be 
valid. 

This paper is organized as follows. In §2, we define the 
coupled spin-pseudospin model that describes the spin 
and charge degrees of freedom of the anisotropic Hub- 
bard ladder at quarter filling. Some details of the method 
of calculation are also given. In §3, we present the results 
of calculation which include the staggered susceptibility 
for pseudospins, the spin and pseudospin excitation spec- 
tra, and the temperature dependence of the uniform spin 
susceptibility. Discussion on the experimental relevance 
to cv'-NaV^Os and summary of the paper will be given 
in 84. 
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2. Model and Method 

Our effective spin-pseudospin Hamiltonian may be 
written as a sum 



H = Ho + Hst 



(1) 



of the quantum Ising Hamiltonian for pseudospins 

^o = Ji(-fE^ + E^ T m) (2) 



and the spin-pseudospin coupling term 

H ST = J 2 ( S * • S *+i - \) C^+i + H - c • ( 3 ) 

i 

The standard notation is used here. Sj and Tj are, re- 
spectively, the spin and pseudospin operators of spin- 1/2 
at site i, where T? = —1/2 (+1/2) means the electron 
is on the left (right) site on the rung of the ladder. J\ 
is the energy scale of the pseudospin system and J2 is 
the coupling strength between the spin and pseudospin 
systems. 

From the second-order perturbation theory, 5, 20, 21 ) we 
have the relations J\ = 2Vj| and J 2 = 4£ 2 /Vj_, where 
t|| and Vj| (t± and V±) are the nearest-neighbor hopping 
parameter and Coulomb repulsion of the leg (rung) of 
the ladder, respectively. We should then have J\ > J2, 
which we assume throughout the present work. We also 
assume the onsite Coulomb repulsion to be U ^00. Rel- 
ative strength of the transverse field to the pseudospins 
is measured by g = 4t±/J\ = 2tj_/Vj| . Note that g in the 
quantum Ising model represents the relative strength of 
the fluctuation of a charge in the rung: if we assume 
one electron in a rung, we have the prefactor gJ\/2 in 
the first term of eq. (2), which is the difference between 
the energies of the bonding and antibonding levels of the 
rung, 2ij_. Thus, if g (or t±) is large the electron is stable 
in the bonding level of the rung, but if g (or t±) is small 
the effect of Vj| easily leads the system to CO. 

We use the conventional world-line QMC method for 
the analysis of the model. We use a 32-site cluster (where 
a site contains a spin and a pseudospin) with periodic 
boundary condition; the cluster-size dependence of the 
calculated results are examined by using clusters up to 
96 sites but we find no significant size dependence in the 
results. Because the model docs not conserve the to- 
tal pseudospin, we have examined a number of ways of 
the spin flips and confirmed that available analytical re- 
sults are reproduced correctly. 22 ) The maximum-entropy 
method is used to calculate the dynamical quantities like 
the spin and pseudospin excitation spectra. 

3. Calculated Results 

3.1 Staggered susceptibility for pseudospins 
The response function is defined as 

XV = Aa ((Sj{-i\)S!) (S3)(S?)) (4) 
Jo 

where Sj(-iX) is the Heisenberg representation of Sj 
and (• • • ) is the canonical average. Xij 1S Fourier trans- 
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Fig. 1. Temperature dependence of the staggered susceptibility for 
pseudospins XtM calculated for the coupled spin-pseudospin 
Hamiltonian. 



formed to the g-dependent susceptibility x{q), which we 
calculate by the QMC method; the q — > limit gives 
the uniform spin susceptibility x(T) and the staggered 
susceptibility is defined as x(o) at q — tt. In the fol- 
lowing, we use the subscripts S and T as in xs(q) and 
XT (9), which stand for the spin and pseudospin degrees 
of freedom, respectively. 

The phase diagram of the quantum Ising model TCq 
is well known; 23 ) at T = there is a long-range order 
for g < 1 (g = 1 is a quantum critical point), which 
corresponds to the zigzag (or 'antifcrromagnetic') CO. 
The calculated staggered susceptibility for pseudospins 
is shown in Fig. 1, where we find that it shows divergent 
behavior at T — > for g < 1. The dispersion relation 
of the pseudospin excitation observed in the calculated 
dynamical structure factor (shown in Fig. 2, see below) 
agrees well with the exact result: 23 ' 



= y"\/l + .9 2 + 2gcosq. 



(5) 



We find in Fig. 1 that the inclusion of the coupling term 
Hst, which introduces the quantum fluctuation via the 
factor T+Tj" , suppresses the divergence. Thus, we may 
say that the inclusion of the spin degrees of freedom in 
the quantum Ising model for pseudospins leads to the 
unstable long-range CO. 

3. 2 Spin and pseudospin excitation spectra 

The dynamical pseudospin structure factor St(q, w) is 
defined as 



S T (q,r) 



1 

N 



E 



-iq(r 2 -r 1 ) 



(T r \(r)T r z 2 (0)) (6) 
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S T (q,r) = - f dujS T (q,uj)K(uj,T) (7) 

JO 

K(lj,t) = e-^ + e-^- T ^ (8) 

where St(<7, t) is the Fourier transform of the imaginary- 
time correlation function. We use the maximum entropy 
method for the inverse Laplace transformation (or an- 
alytical continuation) to obtain Sr(q, w) from St^t). 
The dynamical spin structure factor Ss(q,u>) is similarly 
defined by replacing the pseudospin operator T* with the 
spin operator 

The calculated results for the pseudospin excitation 
spectra at low temperature {k^T = O.IJ2) are shown in 
Fig. 2, where we find that the spectra are under strong 
influence of the spin-pseudospin coupling term J2 . With 
increasing the coupling strength J2/J1 , the peak of the 
pscudospin spectra shifts to higher energies and simul- 
taneously the spectra are broadened. Thus, the lower- 



energy edge of the peak is not affected strongly by the 
coupling strength J 2 , at least when g is large. It seems 
reasonable to suppose that the scattering of the pseu- 
dospin excitations due to spin excitations causes the 
broadening of the spectra. 

The calculated results for the spin excitation spectra at 
low temperature are shown in Fig. 3, where we find that, 
in contrast to the pseudospin spectra, the spin excitation 
spectra change very little; i.e., the peak position, width, 
as well as the shape of the spectra are not affected by 
the parameter Ji when g > 1. When g is small, however, 
the peak position is slightly shifted to lower energies with 
increasing the value of J\ (see Fig. 3 (a)). 

The dispersion relation of the spin and pseudospin ex- 
citations calculated at low temperature are summarized 
in Fig. 4, which are obtained as the momentum depen- 
dence of the peak position of the spectra. For compari- 
son, we show the dispersion of the quantum Ising model 



Tohru NAKAEGAWA and Yukinori Ohta 



9=2.0 



9=2.0 



(a) T=o.U 2 
■ □ : J/J =2.0 
• o :J/J 2 =4.0 

-:quantum Ising model 



pseudospin 




(b) 



T=0.1J 2 
□ :J,/J 2 =2.0 
O :J/J 2 =4.0 
U =0.6x(i/2)sing 




spin 



0.4 



g=1.0 



n/2 
<7 



3 1 



(C) 7=0.1 J 2 

■ □ :J,/J 2 =2.0 
• O : J/J 2=4.0 

'quantum Ising model 



pseudospin 




gooooooooog^ 



g=i.o 



(d) 



1 1 

r=o.u 2 

□ : J/J 2 =2.0 
o :J,/J 2 =4.0 
<o /J =0.6x(n/2)sinq 




n/2 
<7 



7t/2 

9 



Fig. 4. Dispersion relations of the spin (open symbols) and pseu- 
dospin (solid symbols) excitations calculated at k-gT = O.U2. 
Note that the same data at g = 2 (at 3 = 1) are plotted in 
(a) and (b) (in (c) and (d)) in different energy scales Ji and 
J2- The dotted line in (a) and (c) is the dispersion relation for 
the quantum Ising model eq. (5), and that in (b) and (d) is the 
scaled dispersion relation for the ID antiferromagnctic Heisen- 
berg model eq. (9). 



in Fig. 3 (a) and (c); the gap opens when g > 1, which 
is closed at q = it when g — ► 1, leading to the 'antifer- 
romagnctic' long-range order (or zigzag CO). We note 
that the gap remains open irrespective of the value of 
g when we include the coupling term J 2 . In Fig. 4, we 
present the same dispersion relations in a different en- 
ergy scales, i.e., oo q /Ji and w q / 'J 2 . We find that, unless g 
is small, the spin excitation spectra are always inside the 
charge gap, i.e., inside the gap of the pseudospin excita- 
tion spectrum; when the charge gap is large, the energy 
scale of the spin excitations is separated from the high- 
energy charge excitations. With decreasing g, however, 
the energy of the charge excitation decreases at the mo- 
mentum q = 7r to couple with the spin excitations. We 
find in Fig. 4 (b) and (d) that for g > 1 the dispersion of 
the spin excitation spectra scales very well with J 2 ; i.e., 
it does not depend on the value of J\ . The dispersion of 
the calculated spin excitation spectra is fitted well with 
the dispersion of the ID antiferromagnetic Heisenberg 
model 



UqjJi = 0.6 x - sin q 



(9) 



if we include the factor 0.6 as in eq. (9). The factor is 
independent of J\ for g > 1 and at low T. 
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Fig. 5. Temperature dependence of the uniform spin susceptibility 
Xs calculated for the coupled spin-pscudospin Hamiltonian. The 
solid and dotted curves are the uniform susceptibility for the 
system of noninteracting S = 1/2 spins and that for the ID 
antiferromagnctic Heisenberg model, respectively. 



These results suggest that at low temperatures there 
is a parameter region where the spin degrees of freedom 
behaves independently from the pseudospin degrees of 
freedom; it is when g > 1 and the gap of the pseudospin 
excitation spectra is large, inside of which there is a spin 
excitation spectra. Thus, we suggest the validity of the 
decoupling of the coupling term of the Hamiltonian as 

H ST ^J2^ <^ +T m + H - c -) (Si ' Sm - \) (10) 



with 



{T+Tr +1+ K. C )^Q.Q 



(11) 



which leads to the effective Heisenberg-model description 
of the spin degrees of freedom of our model. 

3.3 Uniform spin susceptibility 

To see the validity of the effective Heisenberg-model 
description further, in particular for its temperature de- 
pendence, we calculate the temperature dependence of 
the uniform spin susceptibility for the coupled spin- 
pseudospin Hamiltonian. The results are shown in Fig. 5, 
where comparisons are made with the uniform suscep- 
tibility for the system of free spins and with that for 
the ID antiferromagnetic Heisenberg model. We find 
that the temperature /cbT/J 2 at which J2Xs(T) shows 
a maximum is lower than that of the ID antiferromag- 
netic Heisenberg model; it becomes lower with decreasing 
the value of g or with increasing the value of J\ / J 2 . In 
other words, the deviation from the Heisenberg model is 
large when the quantum fluctuation of the pseudospins 
is small, which occurs when g is small or J\ is large. 

Now, let us analyze the data more precisely. In order 
to do this, we fit the results with the temperature de- 
pendence of the spin susceptibility of the ID antiferro- 
magnetic Heisenberg model, the so-called Bonner-Fisher 
curve; 24 ) i.e., we introduce the T-dependent effective ex- 
change coupling constant J e s(T) and we determine the 
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values so as to fit the calculated uniform spin suscepti- 
bility Xs(T). If the values of J e g thus obtained do not 
depend on T, it follows that the spin degrees freedom of 
our spin-pseudospin model is reduced to a ID Heisenberg 
model 



n, 



spin 



Jeff ^2 ( Si ' Si+1 ~~ l) 



(12) 



at least for the response to the uniform magnetic field. 
The results are shown in Fig. 6. We find that the esti- 
mated value of J e ff(T) is indeed a constant for temper- 
atures below fceT < O.7J2 at g = 2. A crossover tem- 
perature T* (= O.7J2) is thereby defined. The effective 
exchange coupling constant takes a value 



Jeff ~ 0.6J 2 



(13) 



which is consistent with the value estimated from the dis- 
persion relation of the spin excitation spectra (see §3.2). 
We find that also at g = 4 the scaling behavior holds 
up to a higher temperature (fceT iS O.8J2), but with 
the same value of J ff (see Fig. 6 (d)), demonstrating 
the validity of the effective Heisenberg-model descrip- 
tion at T > T* . At g = 1, however, the temperature re- 
gion where J s{T) takes a constant value is already very 
small, although the value is still J c g ~ 0.6J 2 at T ~ 
K, and at g — 0.5, the value of J e ff at T ~ K deviates 
largely from J off = 0.6J 2 (or decreases strongly when 
J1/J2 is large), where the effective Heisenberg- model de- 
scription completely fails. 

We note here that the crossover temperature T* 
roughly scales with J2 rather than Ji, as seen in Fig. 6. 
One might suppose that it should scale with the size of 
the charge gap: i.e., up to temperatures corresponding to 
the energy of the lowest charge excitations, with which 



the pseudospins can excite, the spin excitations may be 
written in terms of the ID antiferromagnetic Heisenberg 
model. However, as we have discussed in §3.2, the size of 
the charge gap shows a rather complicated behavior and 
does not simply scale with either J 2 or J\. The naive 
picture thus does not hold. It may be said however that, 
since there is no other excitations available, the deviation 
from the ID Hciscnbcrg-model description is necessarily 
due to the pseudospin excitations. 

4. Summary and Discussion 

We have calculated the spin and pseudospin excita- 
tion spectra and the temperature dependence of the uni- 
form spin susceptibility of the coupled spin-pseudospin 
Hamiltonian by using the QMC method. We have first 
shown that, when the pseudospin quantum fluctuation 
is large (g > 1), the dispersion relation of the spin exi- 
tation spectra of our model at low temperatures agrees 
well with that of the ID antiferromagnetic Heisenberg 
model with the renormalized effective exchange coupling 
constant J e s = O.6J2 that is independent of the energy 
scale of the pseudospin system J\. Here, the spin exci- 
tation spectra is well inside the charge gap, and thus the 
spin degrees of freedom are separated from the charge 
degrees of freedom. We then have shown that the tem- 
perature dependence of the uniform spin susceptibility 
of our model is well described again by the ID antiferro- 
magnetic Heisenberg model with the same effective ex- 
change coupling constant J e g- = O.6J2. The description 
is valid up to the crossover temperature T* that is re- 
lated to the pseudospin excitations of the system and 
roughly scales with J2 unless the quantum fluctuation of 
the pseudospins is small (g < 1). We have thus demon- 
strated the validity of the effective Heisenberg-model de- 
scription of the coupled spin-pseudospin model for the 
quarter-filled ladders. It then follows that the coupling 
between the spin and pseudospin degrees of freedom, 
which occurs at g < 1, leads to the anomalous spin and 
charge dynamics of the system. 

Although the real material a'-NaV20s is modeled well 
as a 2D trellis-lattice system rather than a ID ladder 
system and thus we need great caution in the direct ap- 
plication of the present results, it may be interesting to 
have a rough idea of the values of the physical parame- 
ters appropriate for a'-NaV205; according to ref., 3 ) we 
have i|| ~ 0.14 eV, t± ~ 0.30 eV, and V\\ ~ Vj_ ~ 0.8 eV, 
which lead to J x ~ 1.6 eV, J 2 ~ 0.10 eV, and g ~ 0.75. 
Thus, the real material may be in the region of g < 1, 
where the spin degrees of freedom are not separated from 
the charge degrees of freedom. The anomalous response 
of the spin degrees of freedom may therefore be expected. 
We would here point out, e.g., that the value of J e ff esti- 
mated from the uniform susceptibility observed in exper- 
iment (which takes the value - 600 - 700 K at T ~ K) 
decreases with increasing temperature, 9 ) which is consis- 
tent with the results of our calculation. The reported 19 ) 
anomalous temperature dependence of the nuclear spin- 
lattice relaxation rate 1 /T\ is also interesting in this re- 
spect. To clarify the dynamics of the spin-charge coupled 
systems near the real CO phase transition, we however 
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need not only to examine the region g < 1 in greater 
detail but also to include the 2D coupling in the present 
model, which we want to leave for future study. 

Because the anomalous charge dynamics has been no- 
ticed also in other transition-metal oxides 25 ) and some 
organic systems, 26 ) we hope that the present study 
will stimulate further researches on the intriguing inter- 
play between the spin and charge degrees of freedom of 
strongly correlated electron systems with CO instability. 
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